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ABSTRACT 

Using solar vector magnetograms of the highest available spatial resolution 
and signal-to-noise ratio we perform a detailed study of electric current patterns 
in two solar active regions: a flaring/eruptive, and a flare-quiet one. We aim to 
determine whether active regions inject non-neutralized (net) electric currents in 
the solar atmosphere, responding to a debate initiated nearly two decades ago 
that remains inconclusive. We find that well-formed, intense magnetic polarity 
inversion lines (PILs) within active regions are the only photospheric magnetic 
structures that support significant net current. More intense PILs seem to im- 
ply stronger non-neutralized current patterns per polarity. This finding revises 
previous works that claim frequent injections of intense non-neutralized currents 
by most active regions appearing in the solar disk but also works that alto- 
gether rule out injection of non-neutralized currents. In agreement with previous 
studies, we also find that magnetically isolated active regions remain globally 
current-balanced. In addition, we confirm and quantify the preference of a given 
magnetic polarity to follow a given sense of electric currents, indicating a domi- 
nant sense of twist in active regions. This coherence effect is more pronounced in 
more compact active regions with stronger PILs and must be of sub-photospheric 
origin. Our results yield a natural explanation of the Lorentz force, invariably 
generating velocity and magnetic shear along strong PILs, thus setting a physical 
context for the observed pre-eruption evolution in solar active regions. 

Subject headings: MHD — Sun: atmosphere — Sun: corona — Sun: flares - 
Sun: surface magnetism — Sun: photosphere 
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Introduction 



In spite of its relative thinness, the solar photosphere plays an important role as the 
transitional layer between the fluid-dominated solar interior and the magnetically dominated 
solar atmosphere. Over the past decades the emergence process of magnetic flux through the 
photosphere has been thoroughly studied by observers and modelers. An early observational 
findi ng was the filamentary, fibril structure of photospheric magnetic flux strands or tubes 



(e.g.. lLivingston fc Harveyill969uHoward fc Stenflolll972tlStenflolll973l ) that defies the spatial 
resolution of our best magnetographs to this day. This fine, fibril stru cture of photospher ic 
magnetic fields is yet to be fully understood and properly interpreted (IParkerl Il979l . 120071 ). 



Soon after the discovery of the filamentary photospheric magnetic fields, it was realized 
that emerg i ng flu x tubes must be twisted in order to maintain their structural integrity 



(Sch uesslerl 1 19 791). Rec ent simulations of magnetic flux emergence ( [Murray et al.l 12006 



Archontis fc Hoodl 120081 and references therein) further showed that magnetic flux tubes 
cannot even emerge without some twist. Twist implies the presence of field- aligned electric 
currents, realized by the curl of the measured or simulated photospheric magnetic field vector 
per Ampere's law. The existence of electric currents automatically implies departure of the 
magnetic structure from its minimum-energy, current-free (potential) configuration, as po- 
tential fields are defined by the gradi ent of a scalar potential calcu lated in a finite, bounded, 
or infinite, partly bounded, volume (jSchmidtlll964j ; ISakurail 119821 ) . Non-potentiality in so- 
lar magnetic fields is the undisputed underlying cause of solar eruptive activity originating 
typically in complex regions of interacting opposite-polarity magnetic fields, interfacing by 
strong polarity inversion lines (PILs). Observed photospheric PILs are typically deformed, 
shredded, and stretched by strong shear flows and associated magnetic shear, measured as 
the local angular differ ence between the observed field vec tor and the expected p otential-field 



config uration (see, e.g.. lWangi (119991 ) forareview, but also lZirin fc Wangl ( 119931 ) . iTiwari et al. 
( 120101 ). and references therein). 



In spite of strong mechanical forces a nd complexity in the ph otosphere, twist seems to 
remain almost unchanged in the corona ( jMcClymont et al.l 119971 ). where the substantially 
decreased plasma /3-parameter allows for a nearly complete spatial filling of the coronal 
volume by magnetic field lines. The low-/? coronal environment further implies almost purely 
field-aligne d currents, except on thin current sheets that interface between flux tubes ((Parker 
19721 . 12004 ). As a result, there is a v olume current in flux tube s that is non-neutralized across 
any cross-section of the corona (e.g.. lLongcope fc Welschll2000l ). The notion of non-neutrality 
of electric currents implies the presence o f a net (nonzero) current for a given m agnetic 
polarity, as defined by previous works (e.g. IWilkinson et al.lll992l ; IWheatlandl 120001 ) . 



On the other hand, the presence of volume currents in twisted magnetic fibrils embedded 
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in a relatively field-free space occupied by fluid plasma, that is the assumed case for the pho- 
tosphere, necessarily implies the existenc e of return (skin or sheath) surface currents at the 
boun daries of these photospheric fibrils (Pizzo 19861 ; Ding et al. 1987 ; Longcope &; Welsch 
2000). In addition to a field-aligned twist component, therefore, photospheric electric current 
density shows a significant component perpendicular to the magnetic field B. This cross- field 
current density partially depends on VB x B-term s, caused by gradients of the magnetic 
field strength B (jZhangil200ll ; iGeorgoulis et al.l 12004 ). As such, this component should peak 
on the surface of the flux tube where the interface between the magnetized and field-free 
mediums occurs. In the simplest case of isolated flux tubes it can be shown (see §5) that 
the Lorentz f or ce due to these cross-field currents is mainly due to magnetic pressure (see, 



e.g., Uacksonl 's (119621 ) decomposition of the Lorentz force into pressure and tension terms) 
and tends to expand the flux tube, hence competing with the twist that keeps it together in 
conjunction with the external non-magnetized plasma pressure. 

Little doubt exists that photospheric magnetic fields are indeed forced i .e., subjected 
to significant Lorentz forces (IMetcalf et al.lll995l ; IGeorgoulis fc LaBontd 12004 ). Above the 
photosphere, however, the plasma pressure decreases abruptly forcing flux tubes to expand 
rapidly to fill nearly the entire coronal volume, thus smoothing out large gradients VB. 
Because of this, return currents cannot reach the corona - they fade when magnetic fields 
become force- free (i.e., showing only a volume current due to the twist along the field). 
Details of this process remain to be determi ned, but an interesti n g me chanism relying on 
torsional Alfven waves has been suggested by lLongcope fc Welsch! (120001 ) . 



While it is generally established that a net electric current exists in an y coron a l cross- 
section, debate pertains over whether photospheric currents are neutralized. iParkerl (j 1996b! ) 
argued that for isolated magnetic fibrils embedded in a field-free volume, regardless of twist, 
the total current must be neutralized across the photospheric cross-section of any given flux 
tube. Simply, then, the vertical compo nents of the volume and the surface currents cancel 
themselves in each fibril. IParkerl (Il996bl ) further suggested that the common practice of infer- 
ring the vertical electric current density J z via the differential form of Ampere's law (V x B) z 
using photospheric magnetic field measurements B bears no physical meaning. Thus, the 
presence of non-vanishing vertical currents in several observational cases is, according to 
Parker, an artifact caused by the limited spatial resolution of the observing magnetograph 
that tends to attribute an artificially continuous magnetic field where magnetic field is in- 
herently discontinuous. 



Melrose (119911 ; Il995l ) followed a different approach, by pointing out that observations 
are more consistent with non-neutralized photospheric current patterns. Therefore, currents 
have to emerge from the solar interior with the emergence of magnetic flux. That V ■ J = 
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implies that currents have to close at the base of the convection zone (i.e., the tachocline) 
where magnetic fields are generated. This has profound implications for flare initiation: 
while Parker's neutralized currents require in-situ storage and release of energy, Melrose's 
large-scale current paths imply the existence of a strong induc tive coupling between the 
coronal flaring volume and the convection zone. iMelrosd (119951 ) described magnetic loops 
as electric circuits (the [E, J] -p aradigm) and cited observations that showed an increase of 
magnetic shear after solar flares (jWang et al.ll 19941 ) to argue that the energy released in flares 
is replenished by pumping non-potential energy through electric currents from deep in the 
solar interior. 



Parkerl (ll996al ) responded to Melrose's proposition by arguing that the (E, J)-paradigm 
results in a set of dynamical equations that are mathematically intractable. In contrast, 
the magnetohydrodynamical (MHD) description of solar processes, in terms of a magnetic 
field B and a velocity field u (the [B, u]-paradigm) is more natural and straightforward. 
The [E, J]-paradigm may hold in laboratory conditions where magnetic fields are generated 
by current-carrying coils wrapped around metals, but not in stellar conditions where the 
dynamo-generated magnetic fields are the ones giving rise to electric currents. 

A number of studies followed or preceded the debate, all focusing on its observa- 
tional aspect, namely, whether there is non- neutralized current in the active-region pho- 
to sphere meaning that ne t electric currents are injected into the atmosphere. The study 
of I Wilkinson et al.l (119921 ) was the only one to conclude that photospheric currents may 
be neutralized. AH other studies (ILeka et al.l Il996l ; ISemel &: Skumanichl Il998t I Wheatland 
2000l ; lFalconerll200ll ) concluded that the active-region photosphere includes current- carrying 
magnetic flux tubes, although currents are roughly balanced at active-region scales, that is, 
non-neutralized currents of opposite senses close onto each other in a given active region. 
The inferred non-neutralized currents are interpreted as sub-photosp heric in origin. W here 
these currents close, however, is unclear from observations and only IWheatlandl (120001 ) ap- 
peared to openly favor Melrose's elect ric circuit anal og. The above studies used both the 
differential form of Ampere's law, that iParkerl (ll996bl ) had criticized for failing to provide a 
physical current, and the integral form of Ampere's law, i.e., 



/ = — I (V x B) • MS = / B • d\ 

4vr Js 4vr Jc 



(1) 



where the total current / is obtained either over a cross-section S of a flux tube, where n is 
the unit vector normal to S, or along a closed curve C bounding the cross-section S. 

Although the expressions in equation (TfJ are equivalent (Stokes' theorem), in practice 
numerical errors can give rise to notable uncertainties. First and foremost, the inference of 
magnetic field components via inversion of Stokes images and profiles includes uncertain- 
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ties subject to the specific inversion model (see, e.g., §1.3 of iBorrero fc Ichimotol (1201 if ) and 
references therein). Because of these uncertainties and since different numbers of pixels are 
being used for the differential and the integral forms of equation (pQ), these uncertainties are 
different for the two theoretically equivalent expressions of Amp'ere's law. The integral rep- 
resentation of Ampere's law is often preferred because it involves fewer magnetogram pixels 
(i.e., the outline of a magnetic polarity, rather than the entire polarity used in the differential 
representation). Moreover, the a zimuthal 180" -ambig;uity, inherent in vector magnetic fields 
measured by the Zeeman effect jHarvevl[l969r ). includes further uncertainties if incorrectly 
resolved. Further , magnetic field measurements are sometimes affected by Faraday rotation 
( ILeka et al.lll996t ISu et al.l 120061 ) that is another source of uncertainties. When evaluating 
equation (JT]), finite differencing introduces truncation errors depending on the various differ- 
entiation schemes used, whereas the curve integral requires precise knowledge of the sequence 
of contiguous points on the bounding contour, regardless of this contour's shape complexity. 

The differential form of Ampere's law gives rise to an electric current density J ~ V x B, 
of which only the vertical component J z can be readily calculated in the photosphere. While 
calculation gives r ise to large J^-val ues, a valid question is whether these current density 
patterns are real. ILeka et al.l ( 119961 ) provide several interesting arguments supporting the 
realism of these currents. The most important is that in case of flux emergence, the increase 
in currents cannot be accounted for by the action of photospheric motions alone. In addition, 
patterns of the vertical current density J z ar e known to reflect th e exp ected mor p holog y 



of the photospheric magnetic field (see, e.g., ISocas-Navarrol ( 120051 ) and iBalthasarl ( 120061 ). 
among others). This evidence, however, cannot be conclus i ve as J z ~ (V x B) z should 
naturally reflect the shape of the field. iMcClymont et al.l (119971 ). in their Appendix B, 
provided further physically-based argum e nts, na mely that (i) the strongest non-neutralized 
currents occur in sunspots, where iParkerl (Il996bl ) admits that the fibril state of the magnetic 
fields probably breaks down, (ii) inclined, spiral flux tubes cannot be maintained in their 
observed position without a true field-aligned current (i.e. twist) reflected in J z , and (iii) 
J z assumes a characteristic pattern along sheared magnetic polarity inversion lines (PILs), 
with narrow ribbon-like structures extending along both sides of the PIL. This, they argued, 
implies that the individual sheared flux strands are indeed blurred due to the insufficient 
instrumental resolution but, nonetheless, the current patterns are real, albeit smoothed and 
onl y partially resolved. An im portant conceptual step along the same lines was achieved 
by ISemel fc Skumanichl (119981 ). who hinted that currents cannot be neutralized in case of 
interacting flux tubes with opposite polarities,that is, when a tight PIL interfaces between 
them. In this case one cannot argue for isolated flux fibrils, an apparent prerequisite for 
Parker's neutralized currents. 



Most of the above studies, including Parker's, underlined the importance of studying 
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vector magnetograms with exceptionally high spatial resolution, su fficient to provide a better 



view o f the fibril, discontinuous magnetic fields. The recent study of lVenkatakrishnan &: Tiwari 



( 120091 ) uses an ultra-high-resolution vector magnetogram to conclude the absence of net cur- 
rents in (isolated) sunspots, as Parker suggested. We also undertake this task here, but 
with a focus on (i) not restricting the analysis to isolated magnetic structures, and (ii) ap- 
plying detailed error propagation to enable safe results and subsequent conclusions. We 
utilize two vector magnetograms obtained by the Spectro-Pol arimeter of the Solar Optical 



Telescope onboard the Hinode spacecraft (IKosugi et al.l 120071 ). As we shall see below, the 



spatial resolution and noise level of the Hinode magnetograms are unparalleled compared to 
magnetograms used in the previously cited studies of the 1990s. Therefore, these data are 
most appropriate and timely for the targeted sensitive calculations. 

The data and preliminary analysis are described in §2. The methodology for calculating 
the total electric currents of "individual" (given the spatial resolution) flux tubes in outlined 
in §3. In §4 we present our results and in §5 we attempt a heuristic physical interpretation 
of them, connecting them with the observed evolution in the photosphere. We summarize 
our study and conclude in §6. 



2. Vector magnetogram data processing 



The Spectro-Polarimeter (SP; iLites et al.l 1200 ll ) is part of the Solar Optical Telescope 
(SOT) package onboard the Hinode spacecraft. The instrument records the four Stokes 
polarization signals at the Fe I 6301.5 and 6302.5 A photospheric magnetically sensitive 
spectral lines by scanning a portion of the solar disk with a spectral sampling of 21.6 mA. 
The normal scanning mode consists of 2048 steps at a nominal pixe l spacing of 0.1585" and 
with a full slit length equal to 1024 pixel lengths ( ILites et al.ll2008l ). 



The inversion of the Stokes profile s is achieved via t he Advanced Stokes Polarimeter 
(ASP) Milne-Eddington inversion code ( ILites et al.l (119931 ) and references therein). The in- 
version returns the magnetic field components with a la sensitivity oBi = 2.4 Mx cm~ 2 (G), 
for the line-of-sight field component, and aB tr = 41 Mx cm~ 2 , for the transverse field compo- 
nent. The full-resolution pixel size being 0.1585 ", tests on quiet- Sun measurements revealed 
a spatial resolution of ~ 0.3" (ILites et al.l 120081 ) . These can hardly be compared with pixel 
sizes of the order 1" or lower (spatial resolution of the order 2" or lower) and noise levels 
aBi, crB tT of the order several tens to 100-200 G, respectively, achieved by magnetographs 
used in the earlier studies of electric current patterns and their neutrality properties. 



The resolution of the azimuthal 180° ambiguity in the selected Hinode vector mag- 
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ne tograms wa s achie ved via the non-potenti al magnetic fie l d calc ulation (NPFC) method 
of iGeorgoulia (120051 ). refin ed as discussed i n iMetcalf et al.l (120061 ). The NPFC method is 
quite insensitive to noise (ILeka et al.l 120091 ) and per forms well also in case of discontinu- 
ous and partially unresolved magnetic flux bundles ( IGeorgoulisI 120121 ) , despite claims by 



Leka et al.l ( 120091 ). The NPFC method is reliable enough to be the method of choice for the 



automatic disambiguation of vector magnetograms obtained by the Vector Spectromagneto- 
graph (VSM) of the Synoptic Optical Long -Term Investigation of the Sun (SOLIS) facility 
rtGeorgoulis et al.ll2008l : Irlennev et al.ll200gh . 



2.1. Hinode/SOT data on NOAA active region 10930 



The first Hinode SOT/SP data set refers to NOAA active region (AR) 10930, a well- 
studied region observed on 2006 December 11. Figure [Jji depicts a continuum image of 
the active- region photosphere, showing a clearly interacting 5-sunspot complex (i.e., inter- 
acting, opposite-polarity pattern including at least one strong PIL). The respective NPFC- 
disambiguated vector magnetogram is shown in Figure [lb. The scanning of the AR started 
at 13:53 UT and ended at 15:15 UT. Even from this single snapshot one notices the extreme 
magnetic shear (~ 90°) along the main PIL of the AR, whose area is selected and expanded 
in the inset of Figure [Tb. Notice also the significant magnetic flux fragmentation everywhere 
in the field of view. Flux strands with cross-section of ~ 1" — 2" are abundant both close 
to the sunspots and to the extended plage areas surrounding them. The typical vertical 
field of these strands ranges between several hundred G and ~ 1.5 kG. Such fragmented 
field concentrations were first observed by the balloon-borne vector magn etograph onboard 



the F lare Genesis Experiment (pixel size 0.18"; spatial resolution ~ 0.5" - iBernasconi et al. 
20021 ). but could only be guessed in previous studies (see, for example, IZirinl (119881 ) for a 
collection of sunspot groups with notable fibril structure observed in Ha wavelengths). Adja- 
cent to the PIL (inset), one notices flux bundles, aligned wit h the PIL with a widt h as small 
as < 2". These "channels" or "lan es" were first reported by IZirin fc Wand (119931 ) and were 



exemplified by IWang et al.l (120081 ) for this particular active region. The total field strength 
of these flux bundles is ~ 2 kG, with a vertical component \B Z \ ~ 1.5 kG. The vertical 
field component on the PIL dropping precipitously to zero, the horizontal field shows an 
amplitude B h ~ (1.5 - 1.9) kG. 

Clearl y, NOAA AR 10930 does not comply with the Hale-Nicholson polarity law (e.g. 
Zirinlll988l ): the leading/trailing polarity pattern is unclear and the AR axis shows a North- 
South, rather than an East- West, orientation. Given the extremely sheared PIL on the other 
hand, it is not surprising that the AR, during its passage from the visible disk (i.e., between 
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12/06/06 and 12/18/06) gave an impressive series of flares (3 X-class, 2 M-class, about a 
dozen C-class, and ~ 100 B-class events). The most intense of these flares, a GOES X3.4, 
took place at 02:14 UT on 2006 December 13, ~ 1.5 days after the Hinode observations of 
Figure [TJ 

A map of the vertical electric current density J z is shown in Figure [lb, with the J z - 
solution for the PIL expanded and shown in the inset. A brief look at the J^-map roughly 
confirms the expected fine structure: in the penumbrae of both sunspots and away from 
the PIL there are roughly r a dial structu r es in agreement with previous studies (for example, 
Semel &: Skumanichl (119981 ): iBalthasarl ( 120061 ). and others). Further, in a secondary PIL 
north of the negative-polarity sunspot there are paired upward/downward J z -patterns along 
the PIL. Along the AR's main PIL (inset) the strong shear indeed creates "tongue"-like 
upward and downward current patterns. Finally, the J^-signature of the extended plage 
fields is much weaker, noisier, and lacks a coherent str ucture, that is evidence of magnetic 
fields nearly normal to the photosphere in these areas ( Martinez Pillet et al.lll997l ). 



Overall, the vertical electric current density solution for NOAA AR 10930 shows strik- 
ing fine structure. Magnetic flux fragmentation easily reaches down to the spatial-resolution 
capability of the observing instrument. Therefore, we cannot claim that the observed mag- 
netic configuratio n is fu lly resolved. The observed /^-structure might, in fact, be what 
McClymont et al.l ( 119971 ) argued for, namely, the result of multiple unresolved magnetic fib- 
rils that are imperfectly observed as a single flux tube. 



2.2. Hinode/SOT data on NOAA active region 10940 

The second data set pertains to Hinode's SOT/SP observations of NOAA AR 10940, 
observed on 2007 February 2. The scanning of the AR started at 01:50 UT and ended at 
02:53 UT. Figure^ shows a continuum photospheric image of the AR. From this image, the 
AR consists of a fairly undisturbed sunspot surrounded by extended plage and a number of 
smaller spots, or pores. Evidence for a small emerging flux sub-region exists just west of the 
main sunspot, where the c ontinuum im age shows dark fibril structures such as those expected 



in arch filament systems ( )Bruzekill967[ ). This area is expanded and shown in the inset. The 
PIL of the AR is in the same area but is much weaker than the main PIL of NOAA AR 
10930. The NPFC-disambiguated solution for the heliographic magnetic field vector is shown 
in Figure [2b. The horizontal field vector shows the typical radial arrangement expected for 
an undisturbed sunspot while some weak shear can be seen along the PIL. 



The solution for the vertical current density, shown in Figure [2b, is also enlightening: the 
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sunspot's penumbra generally exhibits radial J 2 -patterns except close to the PIL, where these 
patterns are deformed (inset). There, the J^-patterns roughly follow the orientation of the 
dark fibrils of the arch filament system as seen in the continuum (Figure [2h). Weak, plage, 
J^-patterns can also be seen, while the small radial pattern in the westernmost part of the 
AR indicates that the darkening in the continuum image and the coherent negative-polarity 
feature in the vector-field image correspond to a small sunspot with almost no umbra. In 
this example, too, the J^-patterns provide physical information on the studied AR. As such, 
they may not be completely artificial. 

The NOAA AR 10940 is much more quiescent compared to NOAA AR 10930 of Figure 
[TJ During its passage from the visible disk (~ 01/28 - 02/09 2007) the AR gave only 3 weak 
C-class flares and ~ 25 B-class events. At the time of the Hinode observations on February 
2, the AR was almost totally flare-quiet. 

We seek the qualitative and quantitative differences between the electric current patterns 
in the two ARs. The findings will be used to determine, using the highest possible spatial 
resolution, whether current-injecting magnetic flux tubes are widespread, rarely occurred, or 
simply non-existent. 

3. The analysis method 

To calculate the total electric current in an entire active region or in polarities of it, 
we take the following steps: first, we translate the pixelized photospheric magnetic-flux 
distribution into a collection of partitions, each with a relatively large 1) number of 
pixels. Second, we apply the integral form of Ampere's law [Equation ([[])] to these partitions 
to calculate their total electric current. Third, we algebraically sum the total currents of all 
partitions to calculate the total current of each polarity and of the entire region. 

Alternatively, the total current can be obtained by, first, deriving the normal component 
J z of the electric current density from a vector magnetogram and then integrating the J z - 
distribution over the area of a given partition with N pixels or the active region as a whole. 
Calculating the total current from the integral Ampere's law does not rely on the current 
densities of all N pixels of a studied partition but it uses the magnetic field components of 
only ~ \/N pixels - those that belong to the perimeter of the partition. Choosing a suffi- 
ciently high magnetic-field threshold in the outline of the partitions we achieve a relatively 
high signal-to-noise ratio, and hence a reasonable uncertainty in our calculations. 
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3.1. Magnetic flux partitioning 



Photospheric flux partitioning is achieved following the partition method of iBarnes et al. 



( 120051 ). Partitioning the vertical component B z of a magnetogram B aims to "discretize" the 
otherwise continuous magnetic flux distribution, by translating it into a series of well-defined, 
non-overlapping partitions of different polarities with uniquely defined outlines, areas, and 
flux contents. Numerically, our method of choice is a "gradient-b ased" flux tessellation 



scheme that uses a downhill-gradient algorithm (e.g. lPress et al.lll992l ) to define photospheric 
partitions based on the two-dimensional morphological and polarity characteristics of B z . In 
addition, a simplification rule is applied, merging partitions that interface by saddle points 
in B z . 

If no thresholds were imposed in the course of partitioning, the entire magnetogram area, 
including active-region and quiet-Sun patches, would be included in an excessive number of 
different partitions. It is of little meaning, however, to include weak and/or tiny parti- 
tions from the highly fragmented Hinode magnetograms, as these partitions are unlikely to 
contribute strong currents. For this reason, we focus on strong-flux active-region patches by 
defining relatively high thresholds for the basic parameters of partitioning: (i) B Zthrea = 100 G 
as the minimum i^-value for each pixel of a given partition, (ii) §thres = 10 19 Mx as the 
minimum magnetic flux $ of a partition, and (iii) S t hres — 40 pixels as the minimum parti- 
tion area S. A partition qualifies for further study only if all three thresholds are exceeded; 
otherwise, it is discarded. 

A slight mean-neighborhood smoothing is also applied to the original i? 2 -maps in order 
to reduce an excessive jaggedness of the partition contours. Once these contours are de- 
termined, however, the original, non-smoothed maps of the horizontal field components are 
used to calculate the total currents. 



3.2. Calculation of the total electric current and uncertainties 

3.2.1. For individual magnetic flux partitions within active regions 

Discretizing the integral Ampere's law [Equation (jl}], we obtain for each qualifying 
partition i the total current 

A Ki 

J ' = 7=E( S -* A ^ + B » A »). ( 2 ) 
fc=i 

where A is the pixel size in physical units and Kj is the number of pixels along the contour 
Ci of the partition. Each k-th pixel of the contour is characterized by a horizontal magnetic 
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field Bh fe = B Xk ic + B Vk y and a vector Al^ = A a . fc x + A^y, whose components assume one of 
three possible values (0, +1, and —1) and represent the unit displacement along the contour. 

A practical issue is the specification of the sequence of contiguous points that com- 
prise the boun ding contour Al o f a pa rtition. For this we use an "edge tracker" algorithm 
introduced by iGeorgoulis et al.l ( 120121 ). The algorithm minimizes the length of the curve 
by iteratively selecting pairs of adjacent c ontour points. C ontour-length minimization is 
achieved by a simulated annealing method ( jPress et al.lll992l ). 



From Equation (T5]), we can estimate the corresponding uncertainty 5Ii of Ii. We assume 
that 5Ii stems only from the uncertainties SB Xk and 5B Vk of the field components B Xk and 
By k , respectively, while A Xk and A yk are known without uncertainties. Considering also 5B Xk 
and SB V to be independent from each other, we obtain from Equation ([2]) 



6L 



cA 

47T 



\ 



E 

k=l 



(AlJBl + Al SB* 



Ilk 



(3) 



If B Xk and B Vk are the horizontal heliographic field components on the image (observer's) 
plane, then their uncertainties, SB Xk and SB Vk , respec tively, can be calculated by im plement- 



ing the analytical method described in Appendix I of IGeorgoulis fc LaBontd (120041 ) . We use 
this approach in the following. 

Thus, for each partition i we calculate the total electric current Ii from Equation (T5]) 
and its uncertainty 5Ii from Equation ([3]). As a sanity check, Equation (T5]) is also used to 
calculate the total electric current I[ using a potential magnetic field vector calculated from 
the disambiguated solution of the photospheric ^-distribution. An exact potential field 
implies that I- = for every i. However, due to various numerical effects, all I- show a 
nonzero value. Comparison between Ii and I- can help assess the significance of the nonzero 
total current Jj of a given partition i and hence the non-neutrality of this partition. 

To calculate the photospheric potential field B p from a given ph otospheri c B z d istribu- 
tion, we apply the analytical Green's function method proposed by ISchmidtl (I1964I ). This 
method is much slower, but quite more accurate - for a given choice of bou ndary conditions - 
than the numerical fast-Fourier-transform method of I Alissandr akid (Il98ll ). In the analytical 



method we assume B p = —Vip, where ip is a smooth scalar potential, i.e., 



Z7T 



(4) 



defined at a given point r = ix + yy of the photospheric plane. B z , on the other hand, is a 
known function of position r' = x'x + y'y (r ^ r'). 
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Taking into account the above uncertainties, we characterize a given partition i as non- 
neutralized if 

\k\ > 5|/;| and >3 5I { . (5) 

In other words, we impose a 5<r significance level toward the total currents I- of the potential 
field and a 3a significance level toward the calculated uncertainties 81 \. 



3.2.2. For active regions as a whole 

In an given active region, let N be the number of qualifying partitions per the criteria 
of §3.1, and Ii be the magnetic flux and total current, respectively, of a given partition 
i. The total currents J + and /_ of the positive and negative polarities, respectively, of the 



region are then given by (e.g. I Wheatlandl 120001 ) 



N 



I + = £ sfh , where sf = |[1 + ^ 



i=l 
N 



(6) 

I- = E s7 J i , where s- = i[l - ^] , 

i=l 



with uncertainties 



N \ V2 

^+ = [E'X . (7) 



,1=1 

V \ V2 



,i=l 



Evidently, significant (i.e., sufficiently larger than applicable uncertainties) J + , J_ imply 
non-neutralized currents for the respective polarities. Moreover, the net current in the AR 
is given by 



N 



E^ = J + +/ - ( 9 ) 



i=l 



and has an uncertainty 



N 



1/2 



51**= \J2 6I n = J 51* +511 ■ (10) 



The solar plasma is considered electrically neutralized, since all the physical processes 
in the Sun are relatively slow. Then, the conservation of the electric charge in the solar 
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atmosphere implies that the electric current density J must be solenoidal, i.e. V ■ J = 0. 
In integral form, this means that the net current through any closed surface, such as the 
photosphere, must be zero. Taking into account that the strongest coronal magnetic fields 
and so the currents are localized within solar ARs, we expect that magnetically isolated 
ARs should be current-balanced (I net ~ 0). In practice, however, there may be large-scale 
magnetic connections between the AR and its surroundings, so the net current may be 
different from zero. This may also happen if the observed area does not include the entire 
AR. 

To jointly characterize imbalances of the total current and the magnetic flux in a given 
AR, we introduce the following dimensionless parameters: the electric-current imbalance and 
the magnetic-flux imbalance, X imb and J-'imb, respectively, where 

x imb = -j^L, (ii) 

EM 

i=l 

|$ + + $_| 

^"imb = "T • (12) 

$ + - $ 

N N 

Here $ + = Yl s t^i anc ^ ®- = S s T^i are the total magnetic fluxes in the positive and 

i=l i=l 

negative polarities of the AR, respectively. We will consider the currents in a given AR to 
be balanced if the inequality 

< ^imb (13) 

is satisfied. This means that the smaller the ratio Timb/J^imb (ideally X imb <C Xi m b <C 1), the 
more confident we are that the AR is current-balanced. 

If the criterion of Equation ( |T3|) is satisfied, it is also meaningful to obtain a uniform, 
normalized degree of neutrality (or non-neutrality, thereof) within each polarity of the AR. 
For this purpose we introduce the global (current) non-neutrality factor 

x ± = -( + ) ri4) 

This dimensionless parameter determines the degree to which different partitions of the same 
polarity show the same sense of total electric current. For strongly non-neutralized, but also 
coherent polarities, whose partitions predominantly show a given sense of non-neutralized 
total current, one obtains X imb <C ~ 1. For very incoherent, or current-neutralized, 
polarities, X im b ~ Z^ n <C 1. The intermediate cases X im b <C X^ < 1 and X im b < X^ < 1 
imply a degree of incoherence and/or non- neutrality for each polarity. 
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4. Current patterns in the selected ARs 

4.1. NOAA AR 10930 

Figure [3^ shows the flux-partitioned vertical magnetic field of NOAA AR 10930. We 
identified a total of N = 531 flux partitions, of which p = 246 show positive polarity 
and n = 285 show negative polarity. To cross-check our calculation, we infer the mean 
vertical electric current density of each partition in two different ways: first, by dividing 
the total current Li of each partition i by the area Si of the partition (J z . = L-JSi) and, 
second, by averaging all the local vertical current densities J Zq within each partition % (J z . = 
(J Zq )), where J Zq have been obtained from the differential expression of Ampere's law. Close 
correlation between these two J 2 -maps implies precise calculation of the bounding contours of 
the partitions. The results from the two mean-J 2 calculations are depicted in Figures and 
[3H, respectively, and are compared in Figure [3b. Evidently there is close agreement between 
the two alternative current density calculations. The generally small discrepancies are mostly 
caused by uncertainties in the calculation of the differential J Zq via finite differencing. This 
action is expected to slightly overestimate the magnitudes | J Zi | of the differential J Zi . Indeed, 
Figure [3b shows {J z ) in tegrai ~ 0.94(J Z ) differential, implying a flatter least-squares best fit 
(black line) than the analytical equality (red line). 

We now focus on the total electric currents of individual partitions. In random order, 
these currents appear in Figure UK We also check which of these currents survive our 
significance test of equations (jSJ). Notably, 120 of the 531 identified partitions satisfy the 
significance criteria. Of them, however, only ~ 25 have strong enough total currents to be 
discernible in Figure (non-neutralized currents shown with red columns and green error 
bars). Regardless, it is evident that there are flux partitions in the AR with non- neutralized 
currents. What is more instructive is the spatial distribution of these partitions, shown in 
the detail of Figure Sb. We find that partitions with the strongest non-neutralized currents 
are adjacent to the PILs of the AR. The two main interacting sunspots (labeled 1 and 2 in 
Figures UK and [lb) have large total currents with values ~ 45 x 10 11 A and ~ —35 x 10 11 A, 
respectively. Partitions adjacent to the secondary PILs of the AR, some of them labeled (3-5), 
have total currents of the order (5 — 10) x 10 11 A. There is virtually no partition with clearly 
non-neutralized total current that is far from the AR's PILs. We conclude, therefore, that 
the intensely flaring NOAA AR 10930 includes strong non-neutralized currents only where 
opposite-polarity flux-tube footprints are close enough to interact, that is, in proximity to 
PILs. 

NOAA AR 10930 as a whole has a net current 7 net = (8 ± 1.6) x 10 11 A. With much 
stronger I + and J_ (±[80 — 90] x 10 11 A), the AR has a total current imbalance X im b — 0.036. 
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This is more than an order of magnitude smaller than the magnetic flux imbalance of the AR, 
J-jmb — 0.4. Consequently, NOAA AR 10930 as a whole can be considered current-balanced. 
This being said, the non-neutrality factor X^ n ~ 0.80, so Z^ n ^> X imb . As a result, besides 
strong non-neutrality in the AR, there is a large degree of coherence in the sense of electric 
current for each polarity: in particular, ~ 80% of partitions of a given polarity show a given 
sense of electric current (positive/negative current for negative/positive magnetic polarities). 
This opposite current-sense preference per polarity achieves an overall current-balanced AR 
so that non-neutralized currents of opposite polarities approximately close onto each other. 

4.2. NOAA AR 10940 

The flux-partitioned vertical magnetic field of NOAA AR 10940 is depicted in Figure [5K 
Here we identified a total of N = 297 partitions, of which p = 156 show positive polarity and 
n = 141 show negative polarity. The average vertical current densities for each partition via 
the integral and the differential calculation are shown in Figures and Eli, respectively, and 
are compared in Figure [5b. For this flare-quiet AR the average current densities are generally 
smaller than those of the flaring NOAA AR 10930. As expected, the differential J z slightly 
overestimates the electric current density. From Figure [5b, (J z )integrai ~ 0.93( J z ) differential 
(black line). 

Important quantitative differences between NOAA ARs 10930 and 10940 are revealed, 
however, when the total currents of partitions are plotted for NOAA AR 10940: Figure [6^ 
provides these currents in random order. The peak partition current here barely exceeds 
3 x 10 11 A that is more than an order of magnitude smaller than the respective peaks in 
NOAA AR 10930 (Figure Hk.). Of the 297 identified partitions, 51 survive the significance 
test of equations ([5]). Of them, only 10 have total currents in excess of 10 11 A. These 
partitions are all close to the weak, loosely formed PIL of the AR (Figure [6b). We conclude 
therefore, that (i) the flare-quiet NOAA AR 10940 also includes non-neutralized electric 
currents, albeit much weaker than those of the flaring NOAA AR 10930, and (ii) these 
weaker, non-neutralized currents are also found exclusively along the AR's weaker PIL. It 
appears that the stronger (more flux-massive) the PIL, the more intense the non-neutralized 
currents it involves. 

NOAA AR 10940 as a whole has a net current 7 net = (3.8 ±1.7) x 10 11 A, comparable 
to that of NOAA AR 10930, although the total currents 7+ and /_ are substantially smaller 
here (±[10 - 15] x 10 11 A). The total current imbalance in the AR is Z imb ~ 0.063, ~ 4.5 
times smaller than the magnetic flux imbalance Jimb — 0.28. Moreover, the non-neutrality 
factor l± n ~ 0.40, so l± n > X imb . This degree of coherence in NOAA AR 10940 is half that 
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of the flaring NOAA AR 10930: different polarities show a significantly weaker preference in 
a given sense of electric currents. 

Summarizing our findings, (i) NOAA AR 10940 as a whole can be characterized current- 
balanced, just like NOAA AR 10930, (ii) NOAA AR 10940 includes much weaker total cur- 
rents than NOAA AR 10930, (iii) both NOAA ARs 10940 and 10930 involve non-neutralized 
current patterns flowing exclusively along their PILs, with the ones for the flare-quiet AR 
10940 (that has a much weaker PIL) being more than an order of magnitude weaker than 
those of the flaring AR 10930, and (iv) NOAA AR 10940 is much less coherent (a factor- 
of-two difference) in terms of preference of electric-current sense per polarity, compared to 
NOAA AR 10930. 



4.3. Dependence of results on partitioning threshold and varying spatial 

resolution 

To obtain the results of §§4.1 and 4.2 we used a fixed, arbitrarily chosen threshold 
B Zthres = 100 G for the magnetic-flux partitioning and we took advantage of the high spatial 
resolution of Hinode's SOT/SP magnetograms. To complete the analysis, we investigate the 
role of both the arbitrary B Zthrea and the spatial resolution in our findings. 

4-3.1. Varying partitioning threshold 

First, we vary the value of B Zthres from 100 G to 1000 G and perform the same analysis 
as in §§4.1 and 4.2. We find that non-neutralized currents in both ARs continue to accumu- 
late in the partitions identified in Figures H] and [61 that is, close to the PILs of the ARs. In 
Figure [7] we show this behavior for NOAA AR 10930. In Figure [Tki, the black curve indicates 
the fraction of the total (unsigned) magnetic flux of the AR that is represented by the parti- 
tioning for each B Zthres -value. Understandably, as more flux is left out of the partitioning for 
increasing B Zthresl this fraction decreases. The magnetic flux imbalance T\ m b and the total 
current imbalance Ii m b for various B Zthres are shown with blue and red curves, respectively. 
For all B Zthres , X imb is always 5-40 times smaller than J-jmb, regardless of the value of B Zthres . 
Therefore, irrespectively of B Zthreg , NOAA AR 10930 includes non-neutralized current pat- 
terns close to its PILs but stays current-balanced globally. The non-neutrality factor 
is substantial at all times (> 0.8) and increases nearly monotonically for increasing B Zthreg . 
Therefore, we can safely conclude that Z imb <C X^ n ~ 1 so the two different polarities in 
NOAA AR 10930 have both strong non-neutralized currents and a clear preference of sense 
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of electric currents, independently of B Zthrea . 

In Figure [TJd the blue curves show 7 + , for the positive-polarity partitions, and 7_ for 
the negative-polarity ones, while the red curve shows the net current J net . That J + < and 
J_ > indicates a predominantly left-handed, counterclockwise twist in the AR. Notice also 
that 1 7+1 and tend to decrease for increasing B Zt , . This indicates that non-neutralized 
currents become weaker if stronger, less interacting, fields are chosen for partitioning. Indeed, 
moving to thresholds B Zthres from several hundred G to 1 kG, we exclude some parts of the 
PILs of the AR. Regardless, the total current balance of the AR does not depend sensitively 
on B Zthres : | /net | is always ~ 3.5 — 4 times smaller than |/ + | and |/_|. 

The same analysis on NOAA AR 10940 is depicted in Figure [8J Because of the much 
weaker currents involved, increasing B Zthrea in this case leads to more striking quantitative 
changes: for B Zthrea ranging between 400 G and 800 G the current imbalance X im b increases to 
~ 0.5 but still remains smaller than the respective flux imbalance Ximb that ranges between 
0.5 and 0.6 for the same B Zthres -range (Figure Ek). For B Zthrea > 800 G, X im b decreases 
significantly to reach ~ 0.2 at B Zthrea = 1000 G. Apparently the combination of weak total 
currents (Figure [8b) and a strongly imbalanced, inhomogeneous magnetic flux distribution 
where the one or the other polarity is favored for different B Zthres -values is the reason for 
these effects. The non- neutrality factor X^ n initially increases and remains high at ~ 0.7 
until B Zthrea ~ 600 G, thereafter decreasing significantly to reach a value ~ 0.4. At all times, 
though, X^ Q > X imb . The polarity currents 7 + and J_ (Figure [8b) peak at ~ 20 x 10 11 A, ~ 5 
times smaller than those in NOAA AR 10930. Their magnitudes also decrease for increasing 
B Zthrea and this causes I net , and hence X im b, to change significantly. Nonetheless, it would 
be biased to claim that whether NOAA AR 10940 is globally current-balanced depends on 
the value of B Zthrea . For reasonable values of this threshold, sufficient to keep a reasonably 
balanced representation of both polarities in the partitioning, the flare-quiet NOAA AR 
10940 has an overall current-balanced structure, includes non-neutralized currents along its 
PIL, and shows some preference in a given sense of current for each polarity, similarly to its 
flaring counterpart NOAA AR 10930. However, this preference is much weaker than in the 
flaring AR, as are the non-neutralized current patterns NOAA AR 10940 involves. 



4-3.2. Varying spatial resolution 

We now study the impact of changing spatial resolution to the current patterns of the 
two ARs. We first bin each magnetogram using different sampling factors and then we 
re-calculate the subsequent current patterns using the analysis of §3 and keeping a fixed 
B Zt , . = 100 G in all cases. For NOAA ARs 10930 and 10940 the results are summarized in 
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Tables [T] and El respectively. 

For both ARs, we reach the following conclusions: 

(1) Decreasing the spatial resolution causes a decrease in the number of the identified 
partitions and a subsequent decrease in the number of partitions that include non- 
neutralized currents, that is, those satisfying the significance criteria of Equations 
§5§. For the coarsest spatial resolution (pixel size ~ 2.5") both NOAA ARs 10930 
and 10940 include only three partitions. The identification of fewer partitions for 
coarser spatial resolution occurs because the morphological properties of the photo- 
spheric magnetograms - where partitioning relies - are suppressed or smoothed out 
when the resolution is reduced. 

(2) Decreasing the spatial resolution leads to a decrease in the total current magnitudes 
|/ + | and |X_|. Even in this case, however, strong non- neutralized current patterns in 
the flaring NOAA AR 10930 are easily discernible. This is not the case for NOAA 
AR 10940, whose current magnitudes fall almost within the uncertainty margins (~ 
10 11 A) for pixel sizes equal to ~ 1.3" and ~ 2.5". That NOAA AR 10930 seems 
to largely maintain its strong non-neutralized currents around its PILs even for very 
coarse resolution justifies previous studies that advocate for non-neutralized currents 



Wheatland 


2000; 


Falconer 


2001) 



Wilkinson et al. 



( 119921 ) who did not find compelling evidence of non-neutralized current patterns in the 



flaring NOAA AR 2372. 

(3) Regarding the ARs' overall current balance, Table [1] suggests that NOAA AR 10930 
remains current-balanced regardless of spatial resolution. For all cases, Xi m b is 5 - 25 
times larger than X im b. NOAA AR 10940 also appears current-balanced except in case 
the pixel size becomes > 1.3". Then, X imb ~ Fimb but with very weak I + and I_. These 
weak, insignificant currents practically invalidate any conclusions about X im b. 

(4) Regarding the non-neutrality factor X^, it appears that binning consistently increases 
the coherence of each polarity in terms of sense of electric currents. This is the case for 
both ARs, although the results are somewhat mired in case of the flare-quiet NOAA 
AR 10940, where X imb ~ X± n for coarser pixel sizes, due to the much weaker - and hence 
more uncertain - currents involved. 



Briefly, we conclude that varying the partitioning threshold and/or altering the spatial 
resolution does not affect our findings of (i) non-neutralized current patterns along PILs, 
and (ii) globally current-balanced ARs. Only very restrictive partitioning or very coarse 
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spatial resolution can affect these findings, and only for flare-quiet ARs lacking tight, well- 
organized PILs. Large partitioning thresholds or severe binning can cause artifacts, such as 
the disa ppearance of these w eak PILs. PILs, however, are parts of the "topological skeleton" 
of ARs (IBungey et al.lll996l ). which means that they are topologically stable features. As 
such, at least the strongest of them will survive, regardless of thresholding or smoothing. 
This study shows that non-neutralized current patterns are exclusively related to PILs; 
apparently the stronger the PIL, the more intense the non-neutralized currents it involves. 
If well-formed PILs cannot disappear due to thresholding or binning, their currents cannot 
disappear under these actions, either. 



5. Physical implications of our results 



We now attempt to establish a qualitative connection between our findings and the 
observed evolution along well-formed PILs. A cruc ial aspect of this evo l ution is magnetic 
shear. Shear is invariably formed along PILs (e.g., iFalconer et al.lll997t lYang et al.l I2004J ) 
with the effect being more pronounced in stronger PILs. Pending confirmation with larger 
active-region samples, this work also implies that non-neutralized electric currents are also 
stronger in stronger PILs. A connection between non-neutralized currents and magnetic 
shear, therefore, appears plausible. 

Another notable finding of this study is that partitions of the same polarity tend to have 
the same sense of non-neutralized currents. This is reflected in the non-neutrality factor Xj^ 
that is ~ 1 for the eruptive NOAA AR 10930, despite the large number (hundreds) of par- 
titi ons for each po l arity. The same result in the same active region was recently reported 
by iRavindra et al.l ( 120111 ). A similar effect, but at a clearly smaller degree, occurs for the 
non-eruptive NOAA AR 10940. This coherence of electric currents, more pronounced in 
more compact ARs with stronger PILs, must be of sub-photospheric origin. Indeed, numer- 



ous simulations of magnetic flux emergence, (e.g.. iTortosa-Andreu fc Moreno- Insertisl 12009 
Archontis fc Hoodll2010l ; ICheung et al.l 120101 ) . all start from a single flux tube with a given 
twist in the convection zone. Upon emergence in the photosphere, the flux tube undergoes 
substantial fragmentation. We show here that in the course of such a fragmentation electric 
currents largely retain the ir sub-photospheric sense. Strong PILs likely indicate a strongly 
twisted, perhaps braided ( jLopez Fuentes et al.ll2003l ). coherent sub-photospheric flux tube, 
while weak PILs suggest a loosely formed sub-photospheric tube with only the necessary 
coherence to survive its emergence in the solar atmosphere. 

But how, and why, does magn etic shear occur along intense PILs? For NOAA AR 
10930, studied here, ISu et al.l ( 120071 ) concluded that the observed shear is due to sunspot 
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rotation and the east- west motion of the emerging positive-polarity sunspot just south of the 
negative-polarity main sunspot of the AR. The main PIL of the AR is between these two 
spots (Figured]). Here we reveal additional information, namely that strong and systematic 
non-neutralized currents are formed along the PIL and only there. These currents suggest 
that the Lorentz force may be the m o st nat ural cause of shear. This is already suggested 
in several works: Manchester fc Low! (120001 ) analytically showed that the tension compo- 
nent of the Lorentz force can cau se shear in a n undu l ated emerging; flux tube. This result 
was numerically demonstrated by iManchesterl (120011 ). Manchester et al.l (120041 ) modeled a 
Lorentz-force-driven shear with a gradient ranging from ph otosphere to corona, su pporting; 
earlier observations of differential shear in active regions (jSchmieder et al.lll996f ). A sys- 
tematic modeling of the Lorentz-force-driven shear that demonstrates a coupling b etwee n 
sub-photospheric and atmospheric magnetic fields was performed by IManchesterl (120071 ). 
This author demonstrated by MHD flux emergence simulations that a gradient in the axial 
magnetic field (along the PIL) during the emergence of an fi-loop is the cause of the Lorentz 
force that further causes shear flows leading to magnetic shear. These flows have an am- 
plitude of about half the local Alfven sp eed, hence i ncreas ing from the photosphere to the 
corona. These results were supported bv lFang et al.l (120101) who further found that Z 7-loop s 
may be formed along PILs, in agreement with iTortosa-Andreu Moreno-Insertisl ( 120091 ). 
In the sub-photospheric part of [/-loops, however, they showed that the Lorentz force can 
unshear the magnetic field lines, contrary to the situation in and above the photosphere. 

In an effort to provide a heuristic explanation of the shear-generating Lorentz force F 
and its physical conn ection to non- neutralized currents, we write F in terms of magnetic 
tension and pressure (|Jacksonlll962l ). i.e., 



±[(B.V)B-Iv B = 



(15) 



The azimuthal component F v of the Lorentz force in cylindrical coordinates (r,(p,z) is given 
by 



4tt dr 



1 d{B% - B 



B 9B « 



+ 



B r B,. 



(16) 



2r dip dz 
All terms in Equation ( lT5"j) are due to magnetic tension, with the exception of the term 
— (l/(2r))dB 2 /d(p that is due to magnetic pressure. 

We first assume a pair of undisturbed, twisted flux-tube footprints embedded in field- 
free, plasma-filled space (Figure IHK)- Obviously there is no strong PIL interfacing between 
the footprints in this case. We further ignore the (generally nonzero) azimuthal Lorentz 
force in the flux-tube interior and focus on the outer edges of the tube (dark blue and red 
rings) where sheath currents reside, and we seek the azimuthal Lorentz force acting in 
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this area. There, the two terms of F v depending on the radial field component B r vanish 
because the field has to be largely azimuthal, with a weaker vertical component B z , to ensure 
a divergence-free magnetic field vector. After some algebra, Equation ffTB]) on the footprint 
edges becomes 

+ (17) 
*- 4vr l r dip + dz } " 10 

Notice that a// derivatives in Equation (TT7|) are due to magnetic tension. Azimuthal sym- 
metry in the undisturbed case implies 8/ dip ~ and since is weak and vanishing at the 
edge of the flux tubes, F v ~ from Equation (TT7j) on the interface between the flux tube and 
the field-free space. On this interface it is the radial component of the Lorentz force that 
becomes important. This component, driven by magnetic pressure, will force the flux tube's 
cross-section to increase as much as the outside plasma pressure allows. The formation of 
roughly azimuthally symmetric - hence fully neutralized within a given pola rity - cross-field 



sheath (return) currents will protect the flux tube from disintegrating. As iParkerl (jl996bl ) 
demonstrated, there are no net currents injected in this situation. 

We now assume a pair of interacting flux-tube footprints that deform due to this in- 
teraction, forming a strong PIL between them (Figure [Hb). The field will still have to be 
predominantly azimuthal on the interface in this case, too, so Equation ( fl71) can still be 
used. However, (i) the vertical field component B z is non- vanishing any more; in fact, it 
can be quite strong along both sides of the PIL, and (ii) enhanced, PIL- localized azimuthal 
field gradients build-up in the deformation area, even in this simplified, axially symmetric 
(with respect to the PIL) case. These gradients preferentially enhance sheath currents along 
the PIL. These currents cannot be neutralized any more within a given polarity. Evidently, 
then, F v ^ along the PIL (but not on it, because B z = there) unless the two terms in 
the parenthesis of Equation (fTTl) cancel each other. This resulting azimuthal Lorentz force 
is purely due to magnetic tension. Cancellation of the derivative terms may conceivably 
happen, but in general it is not the cause, and it is certainly not true in case of asymmetries 
in the magnetic flux distribution within a given footprint. Nonzero azimuthal Lorentz force 
will cause shearing motions along the PIL that will further give rise to magnetic shear, as 
demonstrated in numerical simulations. Qualitatively, this result is similar for all interacting 
opposite-polarity flux tubes. 

In the actual solar photosphere, the plasma /9-parameter is high enough for MHD to 
compete with the hydrodynamics of the non-magnetized plasma. This dense plasma, there- 
fore, should exhibit an inertia to the action of the Lorentz force that can only be overcome 
locally, where strong fields give rise to powerful azimuthal Lorentz forces (Equation (TTTj) ) . 
An intuitive PIL field strength that, if exceeded, should allow F v to move the plasma, is the 
equipartition value B eq needed for /3 = 1. Assuming typical photospheric values for convec- 
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tive flows and mass density, or for number density and the effective temperature, one finds 
B eq ~ (200 - 1400) G, with a mean B eq ~ 800 G. In NOAA AR 10930 the field strength 
around the PIL always exceeds 1.5 kG and, in many cases, it is higher than 2 kG (Figured]). 
In this case it is clear that (3 < 1 in the PIL area, hence F v should be powerful enough to 
move the plasma. In NOAA AR 10940 the PIL field strength ranges between a few hundred 
G to <1.5 kG. In this case, therefore, despite some non-neutralized currents, the azimuthal 
Lorentz force may not be able to shear plasma and magnetic field considerably. This rough 
interpretation qualitatively agrees with the observational facts for both studied ARs. 

In a symmetric case such as the one of Figure [9b, Equation ffTTj) clearly implies an 
opposite orientation of in the two sides of the PIL that leads to an opposite sense of shear 
flows. In case of strong asymmetries, however, the competition of the two tension terms 
in the parenthesis may result in the same sign of F v in the two PIL sides, and hence in a 
similar sense of shear flows. Asymmetries can be simply thought of as severe flux imbalance 
and/or different field morphologies of the partitions deformed along the PIL. For the studied 
NOAA AR 10930, for example, inspection of continuum and magnetogram movies implies 
that velocity shear mainly acts in the counterclockwise direction in both sides of the PIL. 
Figure [Tb shows clearly that the PIL area is strongly asymmetric. Even in case of similar 
shear orientation, however, asymmetries lead to strong velocity gradients across the PIL that 
give rise to strong relative motions in the two sides of the PIL. Albeit in a more complicated 
manner, therefore, the end result of strongly asymmetric shear flows is qualitatively similar 
with that of the symmetric case. The directionality of the shear is determined by the direction 
of that is due to the dominant sense of twist in the flux tube. The consistency of the 
shear-flow orientation for the lifetime of the PIL corroborates our findings that there indeed 
exists a dominant sense of twist in interacting opposite-polarity flux tubes. 

The proposed mechanism of Lorentz-force-driven shear along PILs is, of course, qual- 
itative - careful quantitative studies are nece s sary for all aspects of it to be clarified at 



sufficient detail. The study of iTorok fc Klieml (120031 ) first showed numerically that as soon 



as shearing motions are generated in the PIL area, non-neutralized currents inevitably build 
up in the corona. Shear in that study was generated by placing bipolar flux concentrations 
- twisted by photospheric vortex flows - progressively closer to each o ther: the clo s er the 



concentrations, the larger their non- neutralized total current. Further, ITorok et al.l (120111 ) 



recently analyzed th e evolution of electri c curre nts in an emerging flux-rope simulation. As 



predicted earlier by lLongcope fc Welschl (120001 ) , these authors also found that the electric 



current in the coronal part of the emerging flux rope is, indeed, essentially non-neutralized, 
although the physical mechanism underlying this effect has yet to be clarified. Works like 
the above substantiate our findings for non-neutralized currents in active-region PILs while 
our proposed mechanism, at work in case opposite-polarity photospheric flux concentrations 
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approach close enough to interact, provides a physical context for this effect. 



6. Summary and Conclusions 

Using solar magnetic field measurements of the highest available spatial resolution we 
have studied in detail the electric current patterns of the photospheric magnetic configura- 
tions in two solar active regions; a flaring/eruptive, and a flare-quiet one. Our main ob- 
jective was to determine whether emerging and/or evolving active regions inject significant 
non-neutralized electric currents in the solar atmosphere through the photospheric bound- 
ary. We show that such currents are injected solely and exclusively along the photospheric 
magnetic PILs that accompany the dynamical formation and evolution of active regions. 
In our limited sample of two active regions we find that stronger PILs imply more intense 
non-neutralized currents. This result is robust and insensitive to magnetic-field thresholding 
and the spatial resolution of the studied magnetograms, but has to be confirmed by future 
studies with larger active-region samples. 

Despite current non-neutrality in particular locations within active regions, our study 
confirms that active regions as a whole are current-balanced magnetic structures. This means 
that non-neutralized currents injected on one side of a PIL are roughly counter-balanced by 
non-neutralized currents of the opposite sense that are developed on the other side of the 
PIL. 



An additional finding, in agreement with recent results (IRavindra et al.ll201ll ). is the co- 
herence in the sense of electric currents within a given polarity. We have quantified the effect 
by means of a dimensionless parameter dubbed the global non-neutrality factor (Equation 
( I14p ). We find that this effect is more pronounced in the eruptive active region we study, 
but it is also present in the non-eruptive region (Figures UK; and Tables CD, |2J), despite the 
considerable flux fragmentation in the photosphere. We conclude that this effect is of sub- 
surface origin, which stands in agreement with the results of multiple numerical simulations 
that show substantial fragme ntation of a single b uoyant flux tube in realistic simulations of 



emerging active regions (e.g.. ICheung et al.ll2010l ). 



Our findings seem to put a nearly twenty- year-old debate to rest: amid st claim s ranging 



between no injection of net electric currents in the solar atmosphere, per iParkerl 's (jl996bl ) 
isolated flux-tube picture (see Introduction), and injection of net currents in each and every 
magnetic flux emergence episode, we demonstrate that injection of significant net currents 
does occur but it is as rare - or as frequent - as the appearance of intense PILs in the solar 
photosphere. In this case only, Parker's assumption of isolated flux tubes breaks down and 



-24- 



Melrosel 's (Il99ll Il995l ) proposition of non-neutralized photospheric currents becomes valid. 
It is obviously beyond the scope of this work to argue on the detailed nature of the sub- 
photospheric origin of these currents (i.e., the [E, J]- vs. the [B, u]-paradigm). However, we 
argue that strong non-neutralized currents occur exclusively in case of emerging magnetic 
footprints of opposite polarities that are close enough to interact, i.e. to "feel" each other 
and deform as a result of this interaction. Deformation implies preferential enhancement of 
return currents sheathing these flux tubes at the deformed areas. Because sheath currents are 
perpendicular to the axes of the tubes, these perturbations imply the exertion of a Lorentz 
force that, via magnetic tension, generates shear along the PIL, in terms of both shear flows 
and magnetic shear. The more enhanced the sheath currents developed in more flux-massive, 
tight PILs, the stronger the Lorentz force and the resulting shear, provided that the magnetic 
field strength in the PIL area is greater than the equipartition value (plasma < 1), so that 
the Lorentz force can overcome the hydrodynamic inertia and move the plasma (§5). 

The direction of shear flows is typically opposite in the two sides of the PIL (Equation 
(fT7|) ). The magnetic shear angle reflects the sense of twist of the interacting flux tubes: 
for the studied NOAA AR 10930, negative magnetic polarities associate to positive (non- 
neutralized) currents and vice- versa (Figures [3J @J - this implies a left-handed sense of twist in 
the region. An opposite situation (right-handed twist) in a similar field configuration would 
imply an opposite shear orientation in the region. Notice that the main PIL of NOAA AR 
10930 is also strongly asymmetric, leading to predominant counterclockwise, but similarly 
asymmetric, shear motions on both sides of the PIL. We explain this effect in §5. These 
asymmetric flows roughly lead to the same end result as in the case of oppositely oriented 
shear flows due to significant relative motions in the two sides of the PIL. 

We conclude by emphasizing the need for more detailed qualitative and quantitative 
studies of the causal sequence interacting flux tubes [PIL] - non-neutralized electric currents 
- Lorentz force - shear. Such studies may shed new light into the physics of eruption initi- 
ation and pre-eruption active-region evolution. Besides simulations, the role of the Lorentz 
force appears to become evident even from observations, albeit in post-eruption situations 
( [Hudson et al.ll2008t iFisher et al.ll2012l ). Moreover, the coherence of non-neutralized currents 
per polarity, stemming from a common sub-photospheric origin of these currents, directly 
relates to a dominant sense of twist in eruptive active regions. This, in addition, may trans- 
late to a domina nt sense of magnetic helicity in eruptive active regions, as demonstr ated 
by previous (e.g.. iNindos fc Andrews! I2004J ; lLaBonte et al.ll2007t iGeorgoulis et al.l 120091 ) and 



recent (e.g., iTziotziou et al.ll2012l ) works. Undoubtedly, much remain to be revealed, but 



efforts should focus on putting together as many pieces of the solar eruption puzzle as pos- 
sible in the most meaningful physical interpretation possible. We intend to undertake such 
efforts in future studies. 
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Fig. 1.— Views of NOAA AR 10930, observed by SOT/SP on 2006 December 12, at around 
20:30 UT. (a) Continuum image of the active-region photosphere, (b) The corresponding 
photospheric magnetic field vector on the heliographic plane. A vector length equal to the tic 
mark separation corresponds to a horizontal field strength of 2000 G. (c) The corresponding 
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photospheric vertical electric current density. In all images, a part of the AR, enclosed by 
the blue boxes, has been expanded and shown in the insets. Tic mark separation for the full 
images is 10". Tic mark separation for the insets is 2". A vector length equal to the tic mark 
separation in the inset of b corresponds to a horizontal field strength of 1000 G. North is 
up; west is to the right. 
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Fig. 2.— Same as Figure HJ but for NOAA AR 10940, observed by SOT/SP on 2007 
February 2, around 01:50 UT. A vector length equal to the tic mark separation in the inset 
of b corresponds to a horizontal field strength of 1900 G. North is up; west is to the right. 
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Fig. 3. — Flux and current partitioning of NOAA AR 10930. (a) The partitioned vertical 
magnetic field, saturated at ±2000 G. Different partitions are indicated by their enclosing 
contours, (b) Comparison between the mean vertical electric current densities per partition, 
as calculated using the integral (shown in c) and the differential (shown in d) forms of 
Ampere's law. Both current- density maps (c,d) are saturated at ±15mA/m 2 . The linear 
(Pearson) and rank-order (Spearman) correlation coefficients (cc) between the two current 
density estimates are also given. Tic mark separation in all images is 10". 
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Fig. 4. — Total partition currents in NOAA AR 10930. (a) Total currents for each partition. 
Currents that satisfy the significance criteria of equations (JH]) are shown with red columns 
and green error bars. Five of the largest non-neutralized total currents are encircled and 
labeled (1-5). The total currents J + and J_ for positive- and negative-polarity partitions, 
respectively, the total current imbalance X imb , the magnetic flux imbalance J-i m b, and the 
non-neutrality factor Z^ n are also given, (b) Detail of the active region showing the spatial 
distributions of the partitions that satisfy the criteria of equations (jSJ). Different partitions 
are indicated by different shades of blue and red, according to their total current. Total 
currents are saturated at ±10 12 A. Partitions 1-5, with total currents corresponding to the 
five currents identified in a, are also indicated. Tic mark separation is 10". 
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Fig. 6. — Same as Figure HJ but for NOAA AR 10940. Six partitions with the largest 
non-neutralized total currents are identified here. 



□ □ □ Fraction of magnetic flux represented 




Vertical field threshold Bz thres (G) 

Fig. 7. — Total electric currents in NOAA AR 10930 for different partitioning thresholds 
B Zthres . (a) Curves show the fraction of the magnetic flux represented by the partitioning 
(black), the total current imbalance X im b (red), the magnetic flux imbalance J^mb (blue), and 
the current non-neutrality factor (green) for each B Zthrea . (b) The total currents I + and 
I for the positive and negative polarity, respectively (blue curves) and the net current I net 
(red curve). 
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Fig. 8.— Same as Figure but for NOAA AR 10940. 
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Fig. 9. — Simplified concept for the occurrence of shear-generating Lorentz forces along 
lower-boundary PILs in solar active regions. The situation is qualitatively similar for different 
altitudes, along the normal direction z, until the field becomes space-filling, (a) In case of 
non-interacting flux-tube footprints, hence no well-defined PILs, a lack of azimuthal gradients 
within and on the interface between the flux tube and the field-free space inhibits a shearing 
(azimuthal) Lorentz force F v . (b) In case of interacting, deformed flux tubes, inhomogeneous 
azimuthal gradients are formed (dark shapes), preferentially enhancing sheath currents along 
the PIL (dashed line) that become non-neutralized within a given polarity and give rise to 
an azimuthal, shearing Lorentz force F v in the PIL area. In both panels, the field polarity 
is denoted at the center of each footprint, while red- and blue-color gradients suggest the 
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gradient of B. Spiral curves denote twist in the flux tube, and hence the presence of an 
azimuthal field B v . Each footprint is enclosed by its sheath currents (~ V-B x B - dark red 
and blue rings) for the flux tube to maintain coherence. 
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N 

J, nn 
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J_ 


-^net 


-^imb 


X± 
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J\ imb 


(arcsec) 






(xlO 11 A) 


(xlO 11 A) 


(xlO 11 A) 








0.1585 


531 


120 


-84.9 ± 1.1 


92.9 ± 1.2 


8.1 ± 1.6 


0.036 


0.799 


0.398 


0.317 


449 


76 


-86.1 ±0.8 


89.7 ±0.9 


3.7 ±1.1 


0.018 


0.867 


0.447 


0.634 


65 


19 


-58.7 ±0.6 


70.3 ±0.6 


11.6 ±0.8 


0.086 


0.967 


0.583 


1.268 


12 


6 


-47.9 ± 0.4 


57.2 ±0.6 


9.3 ±0.7 


0.089 


1.000 


0.678 


2.536 


3 


2 


-38.8 ±0.6 


50.9 ±0.6 


12.1 ±0.8 


0.132 


0.978 


0.665 



Table 1: Effect of varying spatial resolution on the total currents of NOAA AR 10930; N is 
the total number of partitions, N nn is the number of non-neutralized partitions. The pixel 
size has been modified by binning the original Hinode SOT/SP magnetogram of the AR. 
Shown are the total currents per polarity J_), the net current (/net); the total-current 
imbalance (Ii m b), the current non- neutrality factor (X^J, and the magnetic- flux imbalance 

(Xi m b). 
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0.317 


297 


45 


14.3 ± 1.5 


-10.46 ± 1 


3.8 ± 1.8 


0.063 


0.403 


0.282 


0.634 
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-7.2 ±0.5 


4.2 ±0.7 


0.156 


0.681 


0.474 


1.268 


10 
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7.0 ±0.6 


-1.2 ±0.4 


5.8 ±0.7 


0.514 


0.640 


0.581 


2.536 


3 


1 


3.9 ±0.5 


2.1 ±0.5 


6.0 ±0.7 


0.580 


0.663 


0.541 



Table 2: Same as Tabled] but for NOAA AR 10940. 



